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Analyses of phenomena exhibiting finite-time decay of quantum entanglement have recently attracted con- 
siderable attention. Such decay is often referred to as sudden vanishing (or sudden death) of entanglement, 
which can be followed by its sudden reappearance (or sudden rebirth). We analyze various finite- time decays 
(for dissipative systems) and analogous periodic vanishings (for unitary systems) of nonclassical correlations 
as described by violations of classical inequalities and the corresponding nonclassicality witnesses (or quan- 
tumness witnesses), which are not necessarily entanglement witnesses. We show that these sudden vanishings 
are universal phenomena and can be observed: (i) not only for two- or multi-mode but also for single-mode 
nonclassical fields, (ii) not solely for dissipative systems, and (iii) at evolution times which are usually different 
from those of sudden vanishings and reappearances of quantum entanglement. 



PACS numbers: 42.50.Xa,03.67.Mn 

I. INTRODUCTION 

Decoherence is a crucial obstacle in practical implementa- 
tions of quantum information processing and quantum state 
engineering. Quantum entanglement is especially fragile to 
decoherence. Yu and Eberly |1] (see also earlier studies in 
Refs. yfl) observed that entanglement decay can occur within 
a finite time. This effect has been referred to as entanglement 
"sudden death" or entanglement sudden vanishing (S V) and it 
can be followed by its sudden reappearance (sudden rebirth — 
SR) J3-01- Reference |Q]] has triggered extensive theoretical 
research on entanglement loss in various systems (for reviews 
see Ref. H). Entanglement sudden vanishing was also exper- 
imentally observed lla-Q]. 

Entanglement SV is often considered to be a new form of 
decay of quantum entanglement, which presumably was not 
previously encountered in the dissipation of other physical 
correlations. Here we would like to point out the general oc- 
currence of sudden finite-time decays and periodic vanishings 
of nonclassical correlations. Namely, the SV and SR effects 
can also be observed during the evolution of entanglement 
witnesses JM1 (for a review see Ref. 11210 and nonclassi- 
cality witnesses (also called quantumness witnesses) Il3l4l9ll 
corresponding to violations of classical inequalities. 

A standard approach to study the SV and SR of quantum 
entanglement is based on the analysis of the time evolution of 
entanglement measures, e.g., the concurrence or, equivalently, 
the negativity or the relative entropy of entanglement 111 211 . For 
a two-qubit system, described by a density matrix p, the con- 
currence C(p) is defined by JH: 

C(p) — max f 0, 2 max X t — XA , (1) 



p(i72 <%) <J2)p*{p2 ® 02) and (72 is the Pauli spin matrix. On 
the other hand, the negativity can be defined as H |2j]] : 

N(p) = max ^0, — 2 min/i^ , (2) 

where pj 's are the eigenvalues of the partial transpose p T and 
factor 2 is chosen for proper scaling, i.e., to get N(p) = 1 for 
Bell's states. 

It is worth noting that not all the SRs and SVs of entan- 
glement and its witnesses can be considered standard: A SR 
should appear only after some finite-evolution time after the 
occurrence of the preceding SV. Specifically, let us now an- 
alyze an example: Both |cosi| and max(0,cosi) vanish at 
7r/2, but only the vanishing of the latter function is associated 
with the proper SV and SR effects. 

Both Eqs. (Q~|) and (0 are given as the maximum of zero 
and some functions, which clearly explains the occurrence of 
SVs if p changes in time. By contrast, SVs do not appear 
for the modified parameters C'(p) = 2 max^ Ai — J^i an d 
N'(p) = — 2minj pj, if Ai and pj have continuous deriva- 
tives in time. 

We deduce that analogous SV and SR effects can be ob- 
served for an arbitrary time-dependent parameter F(t), in 
comparison to some threshold value Fo. From a quantum- 
mechanical point of view, the most interesting parameters 
F are the ones which correspond to classical inequalities 

cl 

F > Fo that can be violated for some nonclassical fields, i.e., 
F < Fq as indicated by the symbol < . On the other hand, 

cl 

the symbol > emphasizes that the corresponding inequality 
must be fulfilled for all classical states. Thus, let us truncate 
such parameter F as follows: 



where the A;'s are the square roots of the eigenvalues of 



F ->■ F = max(0,F - F). 



(3) 
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FIG. 1: (Color online) A simple explanation of how to observe 
the SV and SR of nonclassicality witnesses using, as an exam- 
ple, the unitary evolution of single-mode squeezing in the anhar- 
monic oscillator model given by the Hamiltonian l !58t : (a) normally- 
ordered variances S x<f> (dashed curve) and S op t (solid curve), given 
by Eqs. ( 160b and J6U . (b) truncated normally-ordered variances S X(j> 

(dashed curve) and S'opt (solid curve), given by Eqs. d 1 5 fc and l !19t , 
respectively. Quadrature squeezing occurs if S X4 , < or, equiva- 
lently, if the truncated witness S x , > 0. Principal squeezing occurs 

if S'opt < or if the truncated witness S pt > 0. Here, |ao| 2 = 1/2, 
cj>o = <f> = 0, and So = 0. By including damping, one would ob- 
serve finite-time decays, analogously to the standard sudden decays 
of entanglement. 



A simple illustration of this concept is shown in Fig. 1. For 
brevity, such F and F will be referred to as the untruncated 
and truncated nonclassicality witnesses, respectively. The re- 
definition of the witnesses is a key concept in observing the 
SV and SR effects. 

In the next sections we give general arguments and present 
some specific examples of phenomena and nonclassicality 
witnesses to support our conclusions. 

In Sec. II we recall a definition of nonclassicality and 



present a general method of constructing truncated nonclas- 
sicality witnesses that can exhibit both the S V and SR effects. 
In Sec. Ill, we discuss methods of constructing truncated 
entanglement witnesses. We also give a few simple exam- 
ples of truncated nonclassicality and entanglement witnesses. 
Their evolution in some prototype physical models is studied 
in Sees. IV- VI. We conclude in Sec. VII. 



II. NONCLASSICALITY WITNESSES 

In order to test (and characterize) the nonclassical behav- 
ior of a given state p unambiguously, we use the multimode 
Cahill-Glauber s-parametrized quasiprobability distribution 
(QPD) functions defined for -1 < s < 1 by Ha): 

W^(a) = -Trf p TT^K) ) , (4) 



f^(a k ) = i |exp (a fe r - ajtf + ||£| 2 ) D(0 d 2 £, 

(5) 

D(£) is the displacement operator, a is a complex multivari- 
able (ai, a>2, «m), and M is the number of modes. In spe- 
cial cases (for s = 1,0, —1), the QPD reduces to the stan- 
dard Glauber-Sudarshan P function, Wigner W function, and 
Husimi Q function, respectively. 

A well-known criterion of nonclassicality (or quantumness) 
is based on the P function (see, e.g., Refs. IIZ^l ): 

Definition 1 A state p is considered nonclassical if its 
Glauber-Sudarshan P function is not a classical probability 
density (i.e., it is nonpositive). Otherwise the state p is called 
classical. 

We use this definition of nonclassicality although we are 
aware of its drawbacks (see, e.g., Ref. II24IP . It is also worth 
noting that this definition is often extended by a requirement 
of nonsingularity. That is, a classical P function cannot be 
more singular than Dirac's 5 function. But, in fact, the sin- 
gularity of the P function is implied by its nonpositivity (see, 
e.g., Ref. 0). 

Definition 1 can be equivalently formulated via a complete 
set of nonclassicality witnesses corresponding to violations of 
classical inequalities. Here we apply the method of construct- 
ing nonclassicality witnesses proposed in Refs. lfl3l [3 and 
developed in Refs. 11191 12511 . Alternatively, one can apply an 
approach used by Alicki et al. lfl5l - tT7ll . 

Let us analyze an arbitrary AT-mode operator / = /(a, aJ) 
as a function of the annihilation, a = (0,1,0,2, ...,<zm), and 
creation, a7, operators. The P function enables a direct cal- 
culation of the normally-ordered (denoted by ::) expectation 
values of the Hermitian operator f'f as follows: 

(:/t/ : ) = J d 2 cx\f( a , a *)\ 2 P(cx,cx*). (6) 

Then one can apply another criterion of nonclassicality lfl3l 
H: 
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Criterion 1 A state p is classical if (: f*f :) > for all 
functions f. Conversely, if (: y f :) < for some f then the 
state p is nonclassical. 

These conditions can be compactly written as (: p f :) > 

and (: ftf :) < : 0. By analogy with definitions of entan- 
glement witness (see the following section), the normally- 
ordered Hermitian operator : / '/ : can be referred to as (non- 
linear) nonclassicality (or quantumness) witness [14]. For 
convenience, we call the nonclassicality witness (and also en- 
tanglement witness) not only an observable but also its ex- 
pectation value. Note that the understanding of nonclassi- 
cality witnesses is not strictly limited to operators (see, e.g., 
Refs.OEl). A 

By writing / = J^f c -ih^ where a are arbitrary complex 
numbers, one obtains 



/ f /:) = e« 



fi fj 



(7) 



difference (hi — h 2 ) in two systems. This squeezing occurs if 
the normally-ordered variance 

3 = (: [A(n x - h 2 )} 2 :) (9) 
is negative, where AO = O — (O), with O = hi — fi2- It 

cl 

is a purely nonclassical effect as S > holds for any classical 

cl 

fields. Note that S + So > also holds for any classical fields, 
where So > is a threshold value which can be chosen to be 
arbitrary. Thus, one can analyze a kind of "strong" squeezing 

if S + So < 0. In order to observe the SV and SR of this 
strong squeezing we truncate the squeezing parameter S as 
follows: 



S : 



(0,-(: [A(ni - h 2 )} 2 :) - So) 



> 



0. 



(10) 



By replacing A(hi — n 2 ) by (hi — n-i) in Eq. ©, one can 
consider another normally-ordered witness D 1 resulting from 
the classical inequality 



The normally-ordered moments (: fi fj :) can be grouped 
into the following matrix: 



M ( : n) 



(P) 



flh 
flh 



\(--fVi:) (■■ 



flh 
flh 



fhh 



flfN 



\ 



(8) 



fhfN ■■) j 



We call (nonlinear) nonclassicality witnesses not only 
: y f : and (: f'f :) but also the matrices of normally-ordered 
moments (p) and their functions (e.g., determinants). 
The importance of this approach is motivated by the following 
nonclassicality criterion io[l9ll : 

Criterion 2 A state p is nonclassical if there exists f, such 
that det[Mj n ' ) (p)] is negative. 

Thus, if these nonclassicality witnesses are truncated accord- 
ing to Eq. (O, one can predict infinitely many different kinds 
of SV and SR effects. Note that a given nonclassicality wit- 
ness reveals only some specific and limited properties of non- 
classical states. 

It is worth stressing that nonclassicality witnesses are (usu- 
ally) not measures of nonclassicality. A question arises 
whether S V and SR effects can also be observed for some non- 
classicality measures. Below we give an example of quantum 
dynamics leading to the SV and SR of nonclassicality wit- 
nesses but not of nonclassicality measures. 



A. Examples of truncated nonclassicality witnesses 

To find nontrivial examples of SV and SR of some non- 
classicality witnesses, which are not necessarily entanglement 
witnesses (studied in the following section), we analyze the 
squeezing (or sub-Poisson statistics) of the photon-number 



D' = (: (cihi + c 2 h 2 + c 3 ) 2 :) + |c 4 | 2 > (11) 

assuming real parameters Ck (k — 1,2, 3, 4). In the following, 
we apply 



D = max(0,-<: (hi - h 2 + D Q ) 2 :)) > 0, 



(12) 



which is a special case of D' for (01,02,03,04) — 

(1,-1, A),o). 

So far we have only analyzed two-mode witnesses. Clearly, 
it is also possible to observe the SV and SR during the time 
evolution of multi-mode but also single-mode witnesses of 
nonclassicality. We give only two examples of photon-number 
and quadrature squeezings: 

(i) Single-mode photon-number squeezing (also called sub- 
Poisson photon-number statistics) occurs if Mandel's Q- 

parameter is negative, i.e., (: (An) 2 :)/(: h :) < 0. This 
nonclassical effect can also be described by the truncated wit- 



0,- 



(Ah) 



(13) 



(ii) The standard (So = 0) and strong (S > 0) M-mode 
quadrature squeezing can be defined by 



S x , = (: (Ax^) 2 :) ^ (So), 
or, equivalently, via the truncated squeezing witness 



(14) 



S X4 = max(0, -(: (Ax^) 2 :) - S ) > 0, (15) 

where <f> = (cf>i, 4> 2 , 4>m)- The multimode quadrature op- 
erator is given by B23I1 : 



M 



(16) 



rn—1 



4 



is a sum of single-mode phase-rotated quadratures 

x m (4>m) = a m exp(z0 m ) + a) m exp(-i(j) m ). (17) 

The truncated nonclassicality witness S x ,, given by Eq. ( fT31 >, 
can also be used in a single-mode case. The (^-optimized 
quadrature squeezing is referred to as principal squeezing and 
is defined by the witness Il9ll26ll : 

S op t = minS^ < 0, (18) 



or the truncated witness 



Sopt = max(0, — S'opt - S ) = max5 x , > 0. (19) 

<t> 

Note that all entanglement witnesses are also nonclassicality 
witnesses, but not vice versa. An example of the single-mode 
evolution exhibiting the S V and SR of the nonclassicality wit- 
nesses, corresponding to the quadrature and principal squeez- 
ing, is shown in Fig. 1 for the anharmonic model described in 
Sec. VI. 

Explicit examples of many other two- and multimode non- 
classicality witnesses, corresponding to violations of classical 
inequalities, can be found in, e.g., Refs. Ill flsSjlzl [27l43Qh . 



III. ENTANGLEMENT WITNESSES 

An effective method of constructing entanglement wit- 
nesses can be based on the Shchukin-Vogel entanglement 
criterion l3~lll (or its generalizations i32ll ) for distinguishing 
states with positive partial transposition from those with non- 
positive partial transposition (NPT). 

In analogy to the matrices of normally-ordered moments 
M*p' (p), given by Eq. (|S), one can define the following matrix 
of partially-transposed moments: 



/(/7/i) r (/TA) 1 



(ftfN 
fN 



M f (p T 



(flhi 



(flh) T 



\ 



(20) 



\(fVi) r </iJ 2 ) r ••• (fhf N ) r J 



where / = Ylf c if% f° r arbitrary complex numbers c,, 

(fifj) T = to(fi fjP r ) an d r denotes partial transposition. 
The Shchukin-Vogel entanglement criterion l3ll [3211 can be 
written as: 

Criterion 3 A bipartite state p is NPT if and only if there ex- 
ists f, such that det[Mf(p T )] is negative. 

This criterion resembles Criterion 2 of the nonclassicality. 
Thus, analogously to the nonclassicality witnesses, we refer 
to such matrices Mf(fr) of partially transposed moments 
and their functions (like determinants) as (state-dependent 
nonlinear) entanglement witnesses. It is worth noting that 
according to the original definition, entanglement witnesses 
correspond to observables rather than expectation values [9]: 



An entanglement witness is a Hermitian operator W such 
that tr(W Psep) > for all separable states p sep , while 
tr(W p cn t) < for some entangled states p ent . This con- 
cept was later generalized to nonlinear entanglement wit- 
nesses ITolfTlll . Although our usage of the term entanglement 
witness differs slightly from the original usage, we believe 
that it can improve readability of our paper, while keeping un- 
changed the main idea of entanglement witnesses. 

Here we give only two examples of such entanglement 
witnesses based on Criterion 3. Let us apply the following 
Hillery-Zubairy classical inequalities B33I1 : 



(hih 2 ) > |(ai4)| 2 , (ni){n 2 ) > \(aia 2 )\ 2 , 



(21) 



where hi = a|a, is the photon number operator, and di (aj) is 
the annihilation (creation) operator for mode i = 1,2. Thus, 
we can define the following truncated witnesses 



H = max(0, \{did\)\ 2 - (hih 2 )) *> 0, 



(22) 



H' = max(0, |(aia 2 )| 2 - (ni)(n 2 )) > 0. (23) 

which can be positive only for some entangled states, as 

marked by the symbol °> . These inequalities can be de- 
rived in various ways, e.g., from the Cauchy-Schwarz inequal- 
ity 1331] or from entanglement criteria based on partial transpo- 
sition Oil l32ll . Thus, H and H' are entanglement witnesses, 
so the SV of the concurrence implies also the SV of H(t) and 
H'(t) (if they were nonzero for some evolution times). It is 
worth noting that the inequalities in Eq. (f2TT > are satisfied not 
only by separable states but also by all classical states (marked 

cl 

by > ) since they can be derived from nonclassicality criteria 
based on the P function 1 19]. 

Another simple choice of an entanglement witness can 
be related, e.g., to the violation of Bell's inequality. For 
two-qubit states, a degree of violation of Bell's inequality, 
in its version due to Clauser, Home, Shimony, and Holt 
(CHSH) Q , can be defined as (H H] : 



B\p) 



0, max (u-j 
j<k 



(24) 



where Uj (j = 1, 2, 3) are the eigenvalues of Up — Tj Tp,T p - 
is a real matrix with elements — Tr [p {pi <E) Oj)], and &j 
are Pauli's spin matrices. For brevity, although not precisely, 
B is often referred to as a nonlocality (measure). Analogously 
to the concurrence C, the nonlocality B is defined as the max- 
imum of zero and another quantity, which implies that it is 
possible to observe the SV and SR of B(t) in a dynamical 
scenario. 

If a two-qubit state p violates Bell's inequality then it is 
also entangled, but not vice versa, i.e., there are mixed states 
p (e.g., Werner's states discussed below), for which C(p) > 
and B(p) = 0. Thus, B(p) can be considered as an entan- 
glement witness. The S V of an entanglement measure implies 
the SV of an nonlocality measure (if the latter was nonzero 
at some evolution time). Note that for two-qubit pure states 
B(p) = C(p), so in this case the nonlocality is not only an 
entanglement witness but also an entanglement measure. 
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IV. SUDDEN DECAYS OF NONCLASSICALITY 
WITNESSES FOR NONINTERACTING MODES 

Let us first give a simple example of the environment- 
induced sudden vanishing of the entanglement that is closely 
related to the original idea of finite-time sudden decays. As 
a generalization, we also study sudden vanishings of several 
other nonclassicality witnesses, which occur at times different 
than those for the entanglement vanishing. 

By contrast to the following sections, we analyze the entan- 
glement of two modes (qubits), which are not directly inter- 
acting with each other but only with independent reservoirs. 
Specifically, we describe the SV of the nonclassicality of ini- 
tially entangled states, due to interaction with the reservoirs 
under Markov's approximation, by applying the standard mas- 
ter equation for the reduced density operator p: 

d 



at 



p : 



k=l,2 



n k {2a\p k cik - a k 'a\p - pd k a\ 



(25) 



+{n k + l){2dkpa\ - a\a k p - pa\a k )] - -[Hs,p] 



where -f k are the damping rates, h k are the mean thermal pho- 
ton numbers, n k = {exp[hu> k / (fc^T)] — l} -1 , T is the reser- 
voirs temperature at thermal equilibrium, and k^ is Boltz- 
mann's constant. We assume the reservoirs to be at zero tem- 
perature, so we set n\ — n 2 — 0. The Hamiltonian Hs is 
just the sum of free Hamiltonians for the two noninteracting 
system modes. We solve the master equation by applying the 
Monte Carlo wave function simulation with the collapse op- 



/in k a\ 



[37]. 



erators c lk = ^7(1 + n k )a k and c 2k 

It is worth noting that from the standard physical point of 
view, the quantum entanglement between two systems, and 
the related violation of Bell's inequalities, can be considered if 
the systems are spatially separated and are physically uncou- 
pled 1381 . It is seen that this model (contrary to the models 
studied in the following section) satisfies the second condi- 
tion. 

Our example of the environment-induced sudden vanish- 
ing of quantumness and nonlocality is provided for a system 
coupled to two independent reservoirs. It is worth mentioning 
that common reservoirs in some cases can also enhance entan- 
glement both for two qubits and two modes. This is possible 
due to a mixing mechanism rather than an induced interaction 
among them H39I1 . 

Let us analyze the decoherence of the initial Werner-like 
state defined as [361: 



p m (0) = p|* m )(* 



1 



P 



I. 



(26) 



4 

for < p < 1, m = 1, and = (|00) + |H))/\/2. 

Here, / is the identity operator. Under this initial condition 
the solution of the master equation can be given in the stan- 
dard computational basis as B36I1 : 





hit) 






2p % /5i52 



(+) 









2PV9192 





(27) 



where fcW - {2-g 1 ){2-g 2 )+pg 1 g 2 , h { + ] = ff3 _ fe [2-(l + 
p)g k ], and g k = exp(— ^f k t) for k — 1,2. The concurrence 
and nonlocality decay as follows ll36ll : 

1 



C(t) = max \ 0, -y/gi92[ 2p 



V[2-(l+p)si][2-(l+p)5 2 ])}, (28) 



B 2 (t)=m&x(0,2p 2 gi g 2 -1) 



(29) 



respectively. For comparison, we also calculate the decays 
of the two witnesses of the photon-number-difference correla- 
tions: 



S(t) = max 



0,\(9i+9l + 2P9i92)-S 



D(t) = 



0, 20102(1 +p) 



D 2 - D ( 91 - g 2 ) 



(30) 



(31) 



For simplicity, let us assume now the same reservoir damping 
rate 7, so g\ — g 2 = g. Then, the SV times for the above 
entanglement and nonclassicality witnesses can be different 
from each other as they are given by 



(C) _ 

SV 



I 



■p 



— In , 

7 \2(l-p) 



1 



AS) 



2 7 V 2S 



*sv - 27 ln 



1 



2Dl 



(32) 



(33) 



(34) 



(35) 



The results are shown in Fig. 2 assuming some specific values 
of the damping constant 7 and the initial Werner state p\ (0) 
with parameter p. 

In conclusion, we have given a simple example of the de- 
caying entanglement between two qubits, which are not di- 
rectly interacting with each other, but they are only coupled to 
the environment. We have observed the SVs of the two non- 
classicality witnesses, which are different from the SVs of the 
entanglement and nonlocality measures. 



V. PERIODIC SUDDEN VANISHING OF 
NONCLASSICALITY WITNESSES OF INTERACTING 
MODES 

A. Frequency conversion model 

Here we give an illustrative example of periodic sudden 
vanishing of nonclassicality witnesses during a unitary evo- 
lution of two interacting modes. This is in contrast to the 
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FIG. 2: (Color online) An example of the environment-induced sud- 
den vanishing of the nonclassicality witnesses for two noninteracting 
modes. The damping model is described in Sec. IV for the initial 
Werner-like state pi with p — 0.8. Key: the concurrence C (solid 
curve), nonlocality B (dotted curve), and two witnesses describ- 
ing the photon-number-difference correlations: S (dashed curve) for 
So = 0.03 and D (dot-dashed curve) for D = 0.1. 



states, \tp(0)) = J2ni,n 2 c n 1 .n 2 \ni,n 2 ). The total number of 
photons is a constant of motion, fii(t) + fi2{t) =const. 

An important property of the (undamped) parametric fre- 
quency model is that the nonclassicality of an arbitrary state 
is unchanged during its evolution. By applying the results of 
Refs. lEH-Hl], one can find that the time evolution of the QPD 
for the frequency-converter model, described by Eq. d36b . 
with arbitrary initial fields is simply given by 

W is) ( ai a 2 ,t) = W (s) [^1(01,03, -t), / 8 2 (ai,a 2 ,-t),0], 

(39) 

where 2(0:1, 012, t) are the solutions of the corresponding 
classical equations of motion for the frequency conversion 
model: 

/3i(a\, Oi2,t) = oix cos(k£) — ia 2 sin(/tt), 

(82(0:1, ot2,t) = a 2 cos(Kt) — ia± sm(Kt). (40) 

Equation (|39l means that the two-mode QPD for the model 
discussed is constant along classical trajectories. Thus, if the 
initial fields are nonclassical, their degree of nonclassicality 
(as defined, e.g., in Refs. M44I446I0 remains unchanged at any 
evolution times of the system. But yet we can observe SV and 
SR of entanglement and nonclassicality witnesses as will be 
shown in the following subsections. 



standard analysis of sudden decays applied solely to dissipa- 
tive systems. Note that one can easily include the dissipation 
(as studied, e.g., in the former section) to observe the proper 
finite-time sudden decays and SRs analogous to the standard 
ones. 

As a simple model to study SV and SR, let us study the 
parametric frequency conversion described by the interaction 
Hamiltonian 

H = tiK[a\a 2 exp(— iAutt) + ix\a\ exp(iAwi)], (36) 

which is a prototype Hamiltonian describing two linearly- 
coupled harmonic oscillators. It can be applied to a variety 
of physical phenomena including the process of exchanging 
photons between two optical fields of different frequencies: a 
signal mode with frequency uii and an idler mode with fre- 
quency lj 2 . Then a\ and a 2 are the annihilation operators 
for the signal and idler modes, respectively, and k is the real 
coupling constant. For simplicity, we assume a resonant case 
Au> = uj + uj 2 — uj\ . 

The well-known solutions of the Heisenberg equations of 
motion for the signal, b\ (t), and idler, b 2 (t), modes are given 
by & 

b±(t) = aiCOs(Kt) — ia 2 sin(Kt), 

b 2 {t) = a 2 cos(Kt) — ia\ sin(frf). (37) 

The corresponding solution of the Schrodinger equation is 

W)) = Cni.na 7=f r=f~\ 00 ' (38) 

— vii! vn 2 \ 

711,712 v v 

assuming that the system is initially in a superposition of Fock 



B. Evolution of a pure state 

Let us first analyze the parametric frequency conversion for 
the initial state 1^(0)) = |01). The system evolves, according 
to Eq. ( |3"8i >, into 

\ip(t)) = cos(/ef)|01) -isin(«;t)|10). (41) 

It is a nonclassical state described by the following singular 
(so negative) P function 

f d 

P(a l7 a 2 ,t) = 6[p 1 (a 1 ,a 2 ,t)] 1 + — — - 

V dp 2 (ai,a 2 ,t) 

X ^7^ -)s[f3 2 (a l7 a 2 ,t)}, (42) 

which is given in terms of Dirac's 5 function, its derivative, 
and the solutions of the classical equations of motion, given 
by Eq. (l40l >. Elementary calculations lead to the following 
expressions for the concurrence and nonlocality 

C(t) = B{t) = |sin(2/ci)|, (43) 

the entanglement witness describing the violation of the first 
Hillery-Zubairy inequality 

H{t) = isin 2 (2^), (44) 

and the nonclassicality witnesses for the photon-number- 
difference correlations 

S{t) = max [0,cos 2 (2^) - S ] , (45) 
D(t) = max{0,L» [2cos(2«;i) - D Q }} . (46) 
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FIG. 3: (Color online) Simple examples of the sudden vanishing and 
reappearance of the concurrence and other truncated nonclassicality 
witnesses for two interacting modes. The unitary evolution of the fre- 
quency model is shown assuming: (a) the initial pure state |01) dis- 
cussed in Sec. V B and (b) the initial mixed state, given by Eq. J47t 
with p = 0.8, analyzed in Sec. V C. Key: C (thick solid curve) is 
the concurrence, B (thin solid curve) is the nonlocality; H (dotted 
curve) is the entanglement witness, given by Eq. \22\ , describing the 
violation of the first Hillery-Zubairy inequality; S (dashed curve) for 
So = 1/2 and D (dot-dashed curve) for Do = 1 are nonclassical- 
ity witnesses describing the photon-number-difference correlations, 
which are given by Eqs. UOl and d 1 2b . respectively. From the stan- 
dard point of view, a SR should appear only after some finite evo- 
lution time after the occurrence of the preceding SV. It is seen that 
this condition is satisfied for all the witnesses of the mixed-state evo- 
lution (b), but only for some witnesses of the pure-state evolution 
(a). 



Analogously, one also finds the photon-number sub-Poisson 
statistics of the fields as described by the modified Mandel 
parameters Qi = sin 2 (frf) and Q2 = cos 2 (k£). All these 
nonclassical witnesses exhibit periodic SV and SR effects as 
shown in Fig. 3(a). For example, the S V and SR of the concur- 
rence corresponds to the maximum value of S. Analogously, 
we could observe the out-of-phase SVs and SRs of Mandel's 



Q parameters, which can be clearly understood by recalling 
the classical-like interpretation of two linearly coupled oscil- 
lators when one of them is initially excited (Q2 > 0) and the 
other is unexcited (Qi = 0). During the evolution, the excita- 
tion is transferred periodically between the oscillators. 

One can raise an objection concerning the above example 
that a SV of the concurrence is instantly followed by a SR, so 
they are not the proper SV and SR effects. The same behavior 
is found for the other witnesses including D for Dq = 0, and 
S for So = 0. From the more standard, or more orthodox, 
point of view, a SV (of some witness) should not be instantly 
followed by a SR. By contrast, the SV times differ from the 
SR times for D with Dg > and for S with So > [as shown 
in Fig. 3(a)] that is required in the orthodox approach. 

Other, even more convincing, examples of the SV and SR 
effects can be found by analyzing the evolution of initially 
mixed states as will be shown below. 



C. Evolution of a mixed state 

Let us choose the initial state to be a Werner-like state 
po(0), given by Eq. @ for m = and |* ) = (|01> - 
i|10))/v2. This state evolves as follows 

l-p 



A)(*)=p|*o(*)X*o(*)| 



where 



|* (*)> = -^[/-(t)|01)-i/+(t)|10)] 



(47) 



(48) 



with f±(t) = cos(nt) ± sin(«;i). We find the following evo- 
lutions of the entanglement witnesses and the corresponding 
times of the first SV: 



C{t) = max[0,p|c| - (1 -p)/2] =• / 
B 2 {t) =max[0,p 2 (l + c 2 )-l] => / 



(C) 
sv 



= / 



1-P 
2p 



(49) 



(B) _ 
SV 



= / 



(50) 



H{t) = ^max[0,(pc)Ml-p)] => t£? = / 



V 



(51) 

where f(x) = arccosa;/(2K) and c = cos(2k£). The first SR 
occurs at the time 



(i) 



(52) 



for i = C,B, H. It is seen in Fig. 3(b) for p = 0.8 that the 
first SVs and SRs occur in the following order: 



f (B) 



AH) 

L sv 



AC) 



AB) 
l SR 



AH) 
SR 



AC) 



(53) 



On the other hand, the nonclassicality witnesses D and S, 
given by Eqs. ( fTOb and ( fT2b . respectively, evolve as 

S(t) =max[0,(l-p)/2+p 2 sin 2 (2Ki)-S'o], (54) 
D(t) = max[0, (1 -p)/2 + 2D a psm{2Kt) - L> 2 ]. (55) 
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For Sq — and p < 1, we do not observe a complete van- 
ishing of S(t). For = and p = 1 (which corresponds to 
the initial Bell state), S(t) periodically vanishes to zero and 
instantly increases, so it is not a good example of the SV and 
SR effects. However, for < p < 1 we can observe the 
proper SV and SR effects as shown in Fig. 3(b). The first SVs 
occur at the times 



SV 



AD) _ 



7T / V2&Q+P-1 

tt_ (2Dl+p-l 
An 1 V 



4D p 



(56) 
(57) 



and the first SRs occur at t, 

A D ) 



(S) 
SR 



tt/k — igy and t 



AD) 

SR 



37r/ (2k) — ig V . Note that the first appearances of these wit- 



nesses occur at earlier times, i.e., t = tt/(2k) — i g y for 
i = S,D. It is seen that we can always choose threshold 
values So and Do for any < p < 1 in such a way to observe 
the SVs and SRs of these witnesses for the photon-number- 
difference correlations at arbitrary evolution times also when 
the system is disentangled. 



VI. PERIODIC SUDDEN VANISHING OF 
NONCLASSICALITY WITNESSES FOR A SINGLE MODE 

Finally, let us analyze a single-mode anharmonic oscillator 
described by the interaction Hamiltonian 



H 



1 



(58) 



This is a prototype model of various fundamental phenomena 
including the optical Kerr effect. For simplicity, here we refer 
to this effect only. Under this interaction, the initial coherent 
state |ao) evolves periodically into a nonclassical state 



W)) 



KI 2 /2' 



n=0 



exp 



-n(n — l)r 



In), (59) 



where r is a rescaled time nt. It is worth noting that the Kerr 
state, given by Eq. d59l ), becomes at some evolution times 
a superposition of macroscopically distinguishable two ll47ll 
or more H48TI coherent states, which are often referred to as 
the Schrodinger cat and kitten states, respectively. Among 
many nonclassical intriguing properties of the model (see, 
e.g., Ref. B49I1 and references therein}, the Kerr state exhibits 
high-degree quadrature squeezing (ZoL l50ll . We find that the 
single-mode normally-ordered variance S x<f} of the quadrature 
operator — xi(cf)) = aexp(—icj)) + exp(i0) can be 
compactly written as follows: 

S X4i =2|a | 2 [l + /i2COs(ri 2 +T)-/ 2 i(cosT 2 i + l)] (60) 

in terms of the auxiliary functions defined by tm = 

k\a \ 2 sin(lT)+2((j>~(j>o) and/ fe/ = exp{/c|a | 2 [cos(ir)-l]} 
with ao = \ao\ exp(icf>o). Quadrature squeezing occurs if 

S x . < or, equivalently, if the truncated witness S x , > 0, 



defined by Eq. (T5[ with Eq. ( f60b and <p — <f>. For simplic- 
ity, we set a threshold value So to be zero in this section and 
in Fig. 1. By applying the results of Refs. l26l l50ll . we can 
compactly write the 0-optimized variance S op t describing the 
principal squeezing as follows 

S opt (t) = 2|a | 2 (l -/21 - v722 + /4i- 2/12/21 cost'), 

(61) 

where r' = T12 — T21 + r. Analogously to the former 

squeezing criteria, the principal squeezing occurs if Sopt < 

or if the truncated witness S op t > 0, as given by Eqs. ( fT9] l 
and (16 It . Our results are presented in Fig. 1 for some specific 
amplitude of the initial coherent state. 

Note that the periodic vanishing of the entanglement and 
nonclassicality witnesses, analyzed here and in Sec. V, should 
not be confused with the oscillations of the entanglement mea- 
sures in systems interacting with non-Markovian reservoirs 
(see, e.g., Ref. 115 ill ). The SV and SR effects in such systems 
have different character than studied here. Mazzola et al. ll5"lll 
observed the oscillations in short times, which disappear after 
some finite time and are related to the non-Markovian char- 
acter of the reservoirs. In contrast, in the examples presented 
here, the periodic behavior of the nonclassicality witnesses 
persists as being related to the unitary evolution of the states. 

It is worth stressing again that the aperiodic SV and SR ef- 
fects, which are analogous to the typical sudden decays of the 
entanglement, can be observed by inclusion of the dissipation. 
Assuming Markov's approximation, one can apply the master 
equation, given by Eq. d25l l in a special case for a single mode 
(k = 1). Then the SVs and SRs become aperiodic and the 
final SV occurs after some evolution time, which depends on 
the dissipation. However, the dissipation is not a necessary 
condition for the SV occurrence in this model. 

The SV and SR of the entanglement in two-mode dissi- 
pative coupled Kerr models was studied in Ref. [52]. Here 
we showed that the periodic SV and SR of squeezing can be 
observed even in the single-mode nondissipative Kerr model. 
This example confirms our conclusion of the general occur- 
rence of the SV and SR of nonclassicality witnesses even for 
single-mode undamped systems. 



VII. CONCLUSIONS 

We have applied the concepts of the SV and SR of quantum 
entanglement measures to study the SV and SR of entangle- 
ment and nonclassicality witnesses. 

Our main observations can be summarized as follows: 

(i) SVs can be encountered not only in the dissipation of en- 
tanglement but also of other nonclassical correlation parame- 
ters, related to violations of classical inequalities lfl9l l23ll . 

(ii) SVs occur not only in the dissipation of bipartite or mul- 
tipartite (multimode) interacting or noninteracting systems but 
also in a single-qubit or single-mode systems. Our examples 
include single-mode squeezing of photon number, squeezing 
of quadrature operators B23I1 . and violations of other classical 
inequalities llfM . 
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(iii) Non-dissipative systems, which are initially even in 
pure states, can also exhibit periodic SVs of nonclassical phe- 
nomena and the related nonclassicality witnesses. For in- 
stance, the quadrature squeezing of light in a Kerr medium 
exhibits periodic SVs for some finite periods of time. In order 
to observe the proper finite-time sudden decays analogous to 
the standard sudden decays of entanglement [lj], one should 
add dissipation by coupling such systems to the environment. 
The damping causes irregularity and loss of periodicity of the 
evolution of the nonclassicality witnesses. We can conclude 
that the damping accelerates the occurrence of the first SVs 
but it is not a necessary condition for their occurrence. 

With the help of the nonclassicality criteria fl9l l25ll and 
entanglement criteria 1I3TI l32ll . based on moments of the an- 
nihilation and creation operators, as discussed in Sees. II and 
III, it is possible to construct infinitely many nonclassicality 
and entanglement witnesses. These witnesses, after trunca- 
tion according to Eq. ©, can exhibit the SV and SR effects 
when analyzing their time evolution. 

We hope that these observations might motivate deeper 
analysis of SV and SR of various nonclassicality witnesses 



in specific models and also in experimental scenarios. 



Acknowledgments 

We thank Prof. Ryszard Tanas for discussions. A.M. ac- 
knowledges support from the Polish Ministry of Science 
and Higher Education under Grant No. 2619/B/H03/2010/38. 
X. W. is supported by National Natural Science Foundation of 
China (NSFC) with Grants No. 11025527,No. 10874151, and 
No. 10935010. Y.X.L. was supported by the National Natural 
Science Foundation of China under Grant No. 10975080. F.N. 
acknowledges partial support from the Laboratory of Physi- 
cal Sciences, National Security Agency, Army Research Of- 
fice, DARPA, AFOSR, National Science Foundation Grant 
No. 0726909, JSPS-RFBR Contract No. 09-02-921 14, Grant- 
in-Aid for Scientific Research (S), MEXT Kakenhi on Quan- 
tum Cybernetics, and Funding Program for Innovative R&D 
on S&T (FIRST). 



[1] T. Yu and J. H. Eberly, Phys. Rev. Lett. 93, 140404 (2004). 

[2] K. Zyczkowski, P. Horodecki, M. Horodecki, and R. 
Horodecki, Phys. Rev. A 65, 012101 (2001); A. K. Ra- 
jagopal and R.W. Rendell, Phys. Rev. A 63, 022116 (2001); 
L. Jakobczyk and A. Jamroz, Phys. Lett. A 318, 318 (2003); 
S. Daffer, K. Wodkiewicz, and J. K. Mclver, Phys. Rev. A 67, 
062312 (2003); A. Miranowicz, J. Phys. A 37, 7909 (2004). 

[3] R. Tanas and Z. Ficek, J. Opt. B: Quantum Semiclass. Opt. 6, 
S90 (2004); Z. Ficek and R. Tanas, Phys. Rev. A 74, 024304 

(2006) . 

[4] C.E. Lopez, G. Romero, F. Lastra, E. Solano, and J.C. Retamal, 

Phys. Rev. Lett. 101, 080503 (2008). 
[5] T. Yu and J.H. Eberly, Science 323, 598 (2009); J.H. Eberly, T. 

Yu, Science 316, 555 (2007). 
[6] M.P. Almeida, F. de Melo, M. Hor-Meyll, A. Salles, S.P. Wal- 

born, P.H. Souto Ribeiro, and L. Davidovich, Science 316, 579 

(2007) . 

[7] J. Laurat, K.S. Choi, H. Deng, C.W. Chou, and H.J. Kimble, 

Phys. Rev. Lett. 99, 180504 (2007). 
[8] J.S. Xu, C.F. Li, M. Gong, X.B. Zou, C.H. Shi, G. Chen, and 

G.C. Guo, Phys. Rev. Lett. 104, 100502 (2010). 
[9] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. 
A 223, 1 (1996). 
[10] P. Hyllus and J. Eisert, New J. Phys. 8, 51 (2006). 
[11] O. Giihne and N. Lutkenhaus, Phys. Rev. Lett. 96, 170502 
(2006). 

[12] R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, 

Rev. Mod. Phys. 81, 865 (2009). 
[13] E.V. Shchukin, Th. Richter, and W. Vogel, Phys. Rev. A 71, 

011802(R) (2005). 
[14] J.K. Korbicz, J.I. Cirac, J. Wehr, and M. Lewenstein, Phys. Rev. 

Lett. 94, 153601 (2005). 
[15] R. Alicki and N. Van Ryn, J. Phys. A 41, 062001 (2008); R. 

Alicki, M. Piani and N. Van Ryn, J. Phys. A 41, 495303 (2008). 
[16] M. Zukowski, Phys. Rev. A 79, 024103 (2009). 
[17] S.N. Filippov and VI. Man'ko, Phys. Scr. 79, 055007 (2009). 
[18] A. A. Semenov, DY. Vasylyev, and B.I. Lev, J. Phys. B 39, 905 



(2006). 

[19] A. Miranowicz, M. Bartkowiak, X. Wang, YX. Liu, and F. 

Nori, Phys. Rev. A 82, 013824 (2010). 
[20] W.K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
[21] A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 
[22] K.E. Cahill and R.J. Glauber, Phys. Rev. 177, 1857 (1969). 
[23] W. Vogel and D.G. Welsch, Quantum Optics (Wiley- VCH, 

Weinheim, 2006). 
[24] A. Wunsche, J. Opt. B: Quantum Semiclass. Opt. 6, 159 (2004). 
[25] W. Vogel, Phys. Rev. Lett. 100, 013605 (2008). 
[26] A. Luks, V. Pefinova, and Z. Hradil, Acta Phys. Polon. A74, 

713 (1988). 

[27] M.D. Reid and D.F. Walls, Phys. Rev. A 34, 1260 (1986). 
[28] J. Perina, Quantum Statistics of Linear and Nonlinear Optical 

Phenomena (Reidel, Dortrecht, 1991). 
[29] X. Wang, A. Miranowicz, Y.X. Liu, CP. Sun, and F. Nori, Phys. 

Rev. A 81, 022106 (2010); X. Wang, J. Ma, L. Song, X. Zhang, 

and X. Wang, Phys. Rev. E 82, 056205 (2010). 
[30] A. Miranowicz, J. Bajer, H. Matsueda, M.R.B. Wahiddin, and 

R. Tanas, J. Opt. B: Quantum Semiclass. Opt. 1, 511 (1999); A. 

Miranowicz, H. Matsueda, J. Bajer, M.R.B. Wahiddin, and R. 

Tanas, ibid. 1, 603 (1999). 
[31] E.V. Shchukin and W. Vogel, Phys. Rev. Lett. 95, 230502 

(2005) ; A. Miranowicz and M. Piani, Phys. Rev. Lett. 97, 
058901 (2006). 

[32] A. Miranowicz, M. Piani, P. Horodecki, and R. Horodecki, 

Phys. Rev. A 80, 052303 (2009). 
[33] M. Hillery and M.S. Zubairy, Phys. Rev. Lett. 96, 050503 

(2006) . 

[34] J.F. Clauser, M.A. Home, A. Shimony, and R.A. Holt, Phys. 

Rev. Lett. 23, 880 (1969). 
[35] R. Horodecki, P. Horodecki, and M. Horodecki, Phys. Lett. 

A 200, 340(1995). 
[36] A. Miranowicz, Phys. Lett. A 327, 272 (2004). 
[37] S.M. Tan, J. Opt. B: Quantum Semiclass. Opt. 1, 424 (1999). 
[38] A. Wojcik, Science 301, 1183 (2003). 

[39] F. Benatti, R. Floreanini, and M. Piani, Phys. Rev. Lett. 91, 



10 



070402 (2003); F. Benatti and R. Floreanini, J. Phys. A 39, 2689 
(2006). 

[40] W. Louisell, Radiation and Noise in Quantum Electronics 

(McGraw-Hill, New York, 1964). 
[41] R. J. Glauber, Phys. Lett. 21, 650 (1966). 
[42] L. Mista, Czech. J. Phys. B 19, 443 (1969). 
[43] A. Miranowicz, H. Matsueda, M.R.B. Wahiddin, and J. Bajer, 

Acta Phys. Slov. 48, 293 (1998). 
[44] C.T. Lee, Phys. Rev. A 44, R2775 (1991). 
[45] N. Lutkenhaus and S.M. Barnett, Phys. Rev. A 51, 3340 (1995). 
[46] A. Kenfack and K. Zyczkowski, J. Opt. B: Quantum Semiclass. 

Opt. 6, 396 (2004). 
[47] B. Yurke and D. Staler, Phys. Rev. Lett. 57, 13 (1986). 
[48] A. Miranowicz, R. Tanas, and S. Kielich, Quantum Opt. 2, 253 



(1990) ; R. Tanas, Ts. Gantsog, A. Miranowicz, and S. Kielich, 
J. Opt. Soc. Am. B 8, 1576 (1991). 

[49] R. Tanas, in: Theory of N on- Classical States of Light, eds. V. 
Dodonov and V.I. Man'ko (Taylor and Francis, London, 2003) 
p. 267-312. 

[50] R. Tanas, A. Miranowicz, and S. Kielich, Phys. Rev. A 43, 4014 

(1991) . 

[51] L. Mazzola, S. Maniscalco, J. Piilo, K.A. Suominen, and B.M. 

Garraway, Phys. Rev. A 79, 042302 (2009). 
[52] A. Miranowicz and W. Leonski, J. Phys. B 39, 1683 (2006); A. 

Kowalewska-Kudaszyk and W. Leonski, J. Opt. Soc. Am. B 26, 

1289 (2009); J. Phys. B 43, 205503 (2010). 



